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Abstract 

A novel equality relating the rate of energy dissipation to a degree of violation of 
the fluctuation-response relation (FRR) in nonequilibrium Langevin systems is de- 
scribed. The FRR is a relation between the correlation function of the fluctuations 
and the response function of macroscopic variables. Although it has been established 
that the FRR holds in equilibrium, physical significance of violation of the FRR in 
nonequilibrium systems has been under debate. Recently, the authors have found 
that an extent of the FRR violation is related in a simple equality to the rate of 
energy dissipation into the environment in nonequilibrium Langevin systems. In this 
paper, we fully explain the FRR, the FRR violation, and the new equality with re- 
gard to a Langevin model termed a Brownian motor model, which is considered as a 
simple model of a biological molecular motor. Furthermore, applications of our result 
to experimental studies of molecular motors are discussed, and, as an illustration, we 
predict the value of a new time constant regarding the motion of a KIFlA, which is 
a kind of single-headed kinesin. 

1 Introduction 

A molecular motor is a mechano-chemical enzyme that transduces the free energy derived from 
fuel molecules into the mechanical energy of its motion. Many types of molecular motors are 
known, and they play various roles in living organisms. Here, we overview the variety of known 
molecular motors and manners in which they function in living cells. 

One of the most common molecular motors is an acto-myosin system, which plays the key role 
in muscle contraction jl^ . The acto-myosin system consists of two types of filaments: a myosin 
(thick) filament and an actin (thin) filament. The myosin filament comprises head and tail parts, 
and the former plays an important role in the function of the system. The myosin head has 
enzymatic activity and catalyzes hydrolysis reactions of an adenosine tri-phosphate (ATP) into 
an adenosine di-phosphate (ADP) and an inorganic phosphate (Pi) as ATP + H2O — > ADP + Pi. 
The difference between the free energies of the reactants and products induces a relative sliding 
movement of myosin and actin filaments, and is finally transduced into work performed against 
external load or the viscous resistance of the surrounding medium. In the case of acto-myosin, 
this movement leads to a contraction of a muscle fibril, which is a bundle of acto-myosin filaments. 
Because muscle fibrils are packed in parallel to form a muscle, the contraction of muscle fibrils 
finally causes the contraction of a muscle. In addition to muscle contraction, myosins and actins 
are also known to be involved in migration of motile cells and also in the formation of a contraction 
ring in the case of cell division 

Kinesin is another typical example of a molecular motor 3^. This molecule also hydrolyzes 
ATP molecules and steps along a protein filament, termed a microtubule. One of the main roles 
of this molecule is to carry organelles, such as lipid vesicles, over a long distance in a cell. For 
instance, a kinesin carries an organelle along the axon of a neuron, whose length is sometimes on 
the order of meters. Because it takes an extremely long time for organelles to be transported by 
diffusion alone over such a distance, active transporters, such as kinesins, are necessary to allow 
the tip of an axon communicate with the cell body. 

Although the above two species are linear motors, i.e., the motors that slide on a one- 
dimensional filament, rotating motors have also been found. Fq-Fi ATP synthase is a large 
complex of proteins that is found on an inner membrane of mitochondria This motor is 
known to be composed of two parts: Fq and Fi. Fq is embedded in the membrane and it forms 
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a channel, through which protons (H+) flow according to the difference in their electrochemical 
potential across the membrane. This flow of protons is converted into a rotation of Fo around an 
axis of the molecule. This might be understood by analogy with a windmill. Fi is mechanically 
connected to Fq, and they share the same rotating axis. The rotation of Fq induces that of Fi, 
and the mechanical energy of this rotation is in turn used by Fi to catalyze the synthesis of ATPs 
as ADP + Pi — > ATP + H2O. In this fashion, the electrochemical potential difference of protons 
is transduced to the chemical free energy of ATPs. It is believed that the above scenario can 
be rewound; it is believed that if the concentration of ATPs is much larger than those of their 
hydrolysis products, Fq-Fi is forced to rotate in the opposite direction, and protons are pumped 
against the gradient of their concentrations. Indeed, counter rotation of Fi in the presence of 
surplus ATP has already been observed in experiments jH]. 

Besides these examples, a wide variety of molecular motors are already known, and many 
motors are yet to be discovered. However, they share the same property; they catalyze the 
hydrolysis reaction of ATPs and convert the free energy derived from this reaction into the 
mechanical energy of their motion. Their dimensions are around ten nanometers, and their 
movements are quite stochastic. 

In the past decade, it has become possible to directly observe the stochastic motion of a 
single molecular motor due to developments in experimental techniques [3iniI7|. In the single- 
molecule experiments of this kind, the resolutions in space and time are around one nanometer 
and one millisecond, respectively. These high resolutions enable us to detect every step in a trace 
of a movement of a motor. Using such experimental data, kinetic properties of the stochastic 
movement of molecular motors have been investigated (HI IHl ^3 • 

In contrast to the kinetics, the energetic characteristics of molecular motors remain unclear. 
For example, there is no systematic method at present to determine the amount of free energy 
consumed by a single motor at the single-molecule level. Furthermore, we do not know how 
to determine the amount of the input energy delivered to each mode of the movement of the 
motor, although this information is quite important in order to elucidate the mechanism for the 
motion of the motor. Therefore, it is necessary to establish a systematic method to investigate 
the energetic properties of molecular motors on the basis of single-molecule experiments. 

Before considering how we can achieve this goal, we review the recent theoretical develop- 
ments in modeling a molecular motor. In the past ten years, a number of theoretical models 
for molecular motors have been developed. In particular, a class of nonequilibrium Langevin 
systems has been studied intensively, which are termed Brownian motor models \12\ ll.Sj. Here, 
we describe an often referred example of a Brownian motor model, termed a flashing ratchet, 
in detail [HI E] • This model describes the dynamics of a colloidal particle in one-dimensional 
space. This particle is supposed to have several internal states, and transitions from one state 
to another occur stochastically at certain transition rates. This particle feels a potential energy 
that depends on the internal state of the particle. It has been found in this model that this 
particle will acquire non-zero mean velocity in general, although the particle is subject to the 
periodic potentials alone and no macroscopic potential gradient is present. This is related to the 
fact that the detailed-balance condition is not satisfied when the inter-state transitions occur, 
because the transitions make the probability distribution in a steady state different from the 
canonical (Gibbs-Boltzmann) distribution. When the detailed-balance condition is violated in 
this manner, the system is not in equilibrium, but in a nonequilibrium steady state. For a system 
in a nonequilibrium steady state, a directed motion of the particle is not prohibited by the second 
law of thermodynamics if there is no spatial reflection symmetry in the system. This effect is 
a consequence of Curie's principle and often referred to as the ratchet effect I I . See also the 
discussion in Sec. 121 
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Energetics of Brownian motor models have also been considered. Sekimoto provided a foun- 
dation of energetic considerations of Langevin systems in a general context ^15 . The formulation 
developed by Sekimoto has been used to calculate the energetic efficiency of various Brownian 
motor models 1,1 7! . Parrondo et al. also calculated the rate of input energy for several Brow- 
nian motor models, including the flashing ratchet |18j . According to their method, the rate 
of energy input can be calculated provided that the model parameters such as the profiles of 
periodic potentials and transition rates are given. 

Considering these developments, it seems possible to calculate the rate of energy input of 
a molecular motor, if we adopt a Brownian motor model as a model of the motor molecule 
and determine model parameters from experimental data. Unfortunately, however, it is usually 
difficult to determine all the model parameters on the basis of experimental data because the 
model has too many parameters. Even if we could determine all the model parameters by fitting 
the calculated quantities from the model to experimental data, validity of the values of these 
parameters can hardly be verified. As a consequence, the energetic considerations based on the 
experimentally inferred parameters become unreliable. 

In order to remove this difficulty, the authors recently proposed a method by which one can 
determine the rates of energy input and dissipation of a Brownian motor model without knowing 
the model details ^1 HOI . It has been found that the rate of energy dissipation into the heat bath 
with which the system is in contact is related in an equality to the quantities such as the mean 
velocity, the correlation function of velocity fluctuations and the response function of the velocity 
to an external perturbation to the system. The important point is that this equality involves 
experimentally measurable quantities alone. Therefore, we can determine the dissipation rate 
without referring to the details of the model. Furthermore, since the rate of energy dissipation 
is balanced with that of energy input for a system in a steady state, the rate of energy input 
can also be determined in the same manner. Hence, if the assumption that a Brownian model is 
relevant to a molecular motor is correct, with aid of this result, one may be able to determine the 
rate of energy input for a molecular motor by use of the single-molecule experimental techniques. 

This result is also significant from the viewpoint of nonequilibrium statistical mechanics. 
Our result reveals that the rate of energy dissipation is expressed as an extent of violation of the 
fluctuation-response relation (ERR). The ERR is a relation between the correlation function and 
the response function of a macroscopic variable, which is established for systems in equilibrium 
as the fluctuation- dissipation theorem. The fluctuation-dissipation theorem is proved on the basis 
of a condition termed the detailed-balance condition that is satisfied in equilibrium. In contrast, 
the detailed-balance condition is not satisfied at a nonequilibrium steady state, and this leads to 
violation of the ERR in a system at a nonequilibrium steady state. According to our result, the 
extent of violation of the ERR is directly related to the rate of energy dissipation into the heat 
bath, and it serves as a useful measure to characterize a 'distance' from an equilibrium state. 

In this paper, we describe this result in detail with regard to the hashing ratchet model 
as an example. With respect to the ERR, several observable quantities are introduced for this 
model. Then, it is numerically found that violation of the ERR is related to the energy input and 
dissipation rates in a simple relation. Eurthermore, we mathematically prove these results by 
formulating a path-integral representation for the hashing ratchet model. In addition, we discuss 
applicability of this result to experimental studies of molecular motors. Several technical points 
that would be problematic in experiments are addressed. In order to illustrate our method, 
experimental data for KIElA, which is a kind of kinesins, are analyzed and a parameter related 
to the diffusive motion of the motor is predicted. 

This paper is organized as follows. Eirst, the flashing ratchet model is introduced in Sec. |21 
Then, the fluctuation-dissipation theorem and its violation in this model are described and 
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numerically demonstrated in Sec. El Then, we state our main claim on the relation between 
violation of the FRR and the rate of energy dissipation in Sec. 0] and numerically demonstrate 
it with the flashing ratchet model. In Sec. [3 mathematical proofs of the above results are 
presented. In Sec. El several issues related to experimental applications of our method and 
concluding remarks are presented. 



2 Model 

In this Section, we define the flashing ratchet model (See Fig. Next, several well known 
results on this model will be briefly reviewed. 




Figure 1: Schematic image of the flashing ratchet model. The particle has two internal states, 
and 1, according to which it feels different periodic potentials, Uo{x) and Ui{x), respectively. 
The internal state of the particle is stochastically switched from to 1 and 1 to at transition 
rates ^lio{x) and Qoi{x), respectively. 



Let us consider a Brownian particle suspended in an aqueous solution of temperature T. 
Motion of this particle is restricted on a one-dimensional track and its position at time t is denoted 
by x{t). We suppose this particle is subjected to a force derived from a periodic potential, Ua{x), 
which is switched from one to another according to the internal state of the particle denoted by 
(7{t). In the following, we assume the internal state, cr{t), takes a value either of or 1. a{t) is 
supposed to represent a 'chemical state' of the particle, and the value of a{t) is stochastically 
changed from to 1 and vice versa with probabilities per unit time, i.e., the transition rates, 
Qio(x) and Qoi{x), respectively, which may depend on the current position of the particle. 

The equation of motion for this particle is a Langevin equation with the above ingredients as 

= -7^ - K^,^{x{t))+m + em, (1) 

where m is the mass of the particle and the friction coefficient 7 represents the viscosity effect 
from the environment, ^(t) is Gaussian white noise with mean and variance l^k^T, which 
represents the thermal noise from the environment. As stated above, Ua{x) is an ^-periodic 
potential for o" = 0, 1, and a{t) is a Possion process on {0, 1} with transition rates r2io(x(t)) and 
^Qi{x{t)) for transitions from to 1 and 1 to 0, respectively. The prime denotes differentiation 
with respect to x. The last term in Eq. is added here to represent a perturbation force, which 
is applied from outside, e is a dimensionless parameter that is supposed to be much smaller than 
unity, which assures us of smallness of the external force. 
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As usual, since the inertial effect is negligible for a small particle in a viscous solvent, we 
consider the overdamped version of Eq. in the following: 



7i(t) = -f/;(t)(x(t)) + ew+e/pW, 



(2) 



where the dot denotes the time derivative: x{t) = dx{t)/dt. 

For the model given in Eq. (jSJ, by use of standard techniques ^21,, it can be proved that the 
probability distribution Pf^(x,t), which represents the probability to find the system in the state 
a and the position x at time t, obeys the following Fokker-Planck equation |9l I14j. 



d_ 
dt 



Po{x,t) 
Piix,t) 



Qoi{x) Ci{t)-nio{x) 



Po{x,t) 
Pi{x,t) 



where 



ld_ 

7 dx 



UUx)-efpit) + kBT 



d_ 

dx 



(3) 



(4) 



is the Fokker-Planck operator for each state. 

The steady-state distribution in the absence of perturbation, P^^{x), is obtained from Eq. Q 
as a stationary solution with e = 0. It has been found that the canonical distribution 



Kix) = ^exp( 



-U^ix)/kBT), 



(5) 



where Z = Jq Y^l.=Qew{~Ua{x)/kBT)dx, provides the steady-state distribution if and only if 
the transition rates satisfy the condition 



Qoiix) 



-[Uiix)-Uoix)]/ksT 



(6) 



In this case, as explained in Sec. 15.2. ![ the detailed-balance condition is satisfied, and the system 
is in equilibrium. In general, the steady-state distribution differs from the canonical distribution 
when Eq. © is not satisfied, and the system is at a nonequilibrium steady state. In this case, 
the detailed-balance condition is violated, which results in the non-vanishing velocity. 
Analysis of Eq. Q reveals that the mean velocity with e = is expressed as |14) 



1 

7 Jo 



[Pt{x)U;,{x) + Pl\x)U[{x)] dx. 



(7) 



From this expression, it is easily confirmed that the steady velocity v vanishes in the case of 
equilibrium, where equals the canonical distribution P^{x). In nonequilibrium steady 

states, since P^i^x^ is different from the canonical distribution, v differs from zero in general. 

Next, let us consider energetics of this model. In the following discussion, we set e = for 
a while. As long as the particle is moving on a single potential, i.e., it stays in a single state, 
no net energy supply from outside is required because the particle is subject to the potential 
(conservative) force alone. Injection of energy is necessary when the particle undergoes a tran- 
sition from one state a to another state a'. In this occasion, the difference in potential energies, 
Ua-'ix{t)) — Ua{x{t)), must be externally supplied. Thus, by representing the timings at which 
transitions occur during time interval from to t as {ti,t2 ■ ■ ■ ,tn}, the external energy input 
during this interval is represented as 



n 

E{t) = Y,[u.(t,+oMt,))-U, 



{x{t,)) 



(8) 
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which is a stochastic variable. The average energy input per unit time is calculated as 



J . lim M 

i^oo t 

= / [Ui{x)-Uo{x)][ni^{x)P^\x)-n^i{x)Pl\x)]dx. (9) 

The second line follows from the consideration that the probability per unit time at which a 
transition from a to a' at x occurs is equal to Q(yi(j{x)Pp{x) at a steady state. As mentioned in 
Sec.n calculation of input and dissipated energy, J, using Eq. Q requires the detailed knowledge 
on the potential profiles and the transition rates. 

The above definition of input energy is compatible with Sekimoto's argument on energetics 
of Langevin systems. Sekimoto suggested that dissipated energy into the environment during 
time interval from to t can be expressed as 

Q{t)= f [^x{s)-i{s)]odx{s), (10) 

JO 

where the symbol o means that the above Stieltjes integral should be interpreted according to 
Stratonovich's prescription [22]. Q{t) is also a stochastic variable. With aid of Eqs. (|2j and 
we obtain the following relation 

E{t) = Q{t) + {x{t)) - C/,(t) (x(0)), (11) 

which expresses the conservation law of energy. Since \U„{t){x{t)) — [/^-(o) (x(0))] /t ^ holds as 
t — > oo at a steady state, the average rate of energy input is equal to that of energy dissipation 
as limt^oo E{t)/t = limt^oo Q{t)/t = J. 

3 The Fluctuation-Response Relation 

In this Section, we introduce the fluctuation-response relation (FRR) with regard to the flashing 
ratchet model. First, after defining several quantities, the fluctuation-dissipation theorem is 
stated. We then demonstrate the FRR numerically for the flashing ratchet model in equilibrium. 
Next, we numerically show that the FRR is violated in a nonequilibrium steady state of the 
flashing ratchet model. 

3.1 Equilibrium case: The Fluctuation-Dissipation Theorem 

First, let us define several quantities that are involved in the statement of the theorem. One of 
the important quantities is the velocity correlation function defined as follows. 

C{t)^{[±{t)-v][x{Q)-v])^, (12) 

where {■ ■ ■)q denotes the ensemble average in the absence of a perturbation (e = 0). According 
to Wiener-Khinchin theorem, the Fourier transform of C (t) defined as 

/oo 
C(t)e*^Mt, (13) 
-oo 

is equal to the power spectrum density of the velocity fluctuations, x{t) — v. Because a system 
in a steady state possesses the time translational invariance, we find a symmetry property as 
C{t) = C{—t). The imaginary part of C{lo) vanishes because of this symmetry. 
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Next, we define the response function that characterizes the hnear response of the velocity 
to a perturbation force as 

lim = r R^t - s)fp{s)ds, (14) 

where (• • denotes the ensemble average taken in the presence of the perturbation. By adopting 
the perturbation fp{s) as an impulse at time 0, i.e., fp{s) = 5{s), it is found that the response 
function R{t) describes the response of the mean velocity to this impulse as 

ii(t) = lim ^^^^^^^-^ with fp{t) = 5{t). (15) 

Because of the causality, R{t) for negative t is always zero. 

Due to the convolution theorem, the Fourier transform of the response function, termed 
susceptibility, 

poo 

R{uj) = / i?(t)e''^*dt (16) 







1 rim.-'^^r.t. 



satisfies the following relation 

R{uj) = -J— hm / ^Z^I^^^e-Mt, (17) 

where fp{uj) is the Fourier transform of fp{t). This relation enables us to determine the suscep- 
tibility experimentally by applying a sinusoidal perturbation: fp{t) = sinojt. 

Now, we can state the fiuctuation-dissipation theorem. Here, we suppose that Eq. © IS 
satisfied and thus the system without perturbation (e = 0) is in equilibrium. In this case, as 
mentioned above, the mean velocity of the particle vanishes: v = 0. Furthermore, the correlation 
function and the response function satisfy the following relation 

C{t) = kBTR{t) for t > 0, (18) 

or equivalently 

C{uj) = 2kBTR+{uj), (19) 

where i?^(w) represents the real part of R{uj). Proof of these relations is provided in Sec. El In 
what follows, we term Eq. H19|) the fluctuation-response relation (ERR). 

It should be noted that the ERR exactly holds independent of the model datails such as the 
profiles of periodic potentials and the transition rates. As explained in Sec. 15.21 the ERR is 
a direct consequence of the detailed-balance property of a system in equilibrium. Due to this 
simple relation, the measurements of the correlation function and the response function give the 
identical information on the system if it is in equilibrium. 

Now, we numerically demonstrate the ERR for the flashing ratchet model. For the numerical 
calculation, we adopted the model parameters as Uq = sin(27r2;/^) + 1/2 sin(47rx/^) and Ui = 
const., Qio{x) = Qoi{x) = ^. All the quantities are made dimensionless by normalizing 
£ = 7 = k^T = 1. On the basis of model Eqs. @ and ©, we calculated the power spectrum 
density C{u)) and the susceptibility R{uj) by using the Fourier-modes expansion [2211^. In Eig.|2l 
(a), C{ui) and 2kBTR~^{Lo) are plotted as functions of the frequency to. As clearly seen, the ERR 
holds irrespective of the frequency, implying Eq. (|19j) . Since Eq. (|18j) can be obtained by inverse 
Fourier transform of Eq. H19|l . Eq. H18|) has also been verifled. 

^This is a special case of Eq. These transition rates are obtained by choosing ftio(x) = f2e'"^°'^^ and 

noiix) = J^e*^^!'^' and letting Q ^ 0. 
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Figure 2: C{u:) and 2kBTR'^{uj) as functions of lo. (a) Equilibrium case (Qio(x) = Q.qi{x) = 0). 
(b) Nonequilibrium case (J7io(a;) = Q.m{x) = 10). The solid and dotted curves represent C{ijj) 
and 2k^TR'^{uj), respectively. The arrows donote the decay rate uj* . All the quantities are 
dimensionless under the normalization as ^ = ^ = k-^T = 1. 
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We can obtain several important quantities from this plot. First, the diffusion coefficient, 
defined as 

i— >oo 2t 

satisfies D = C(0)/2. Therefore, we can read the value of D from the low-frequency plateau of 
C(u;). Next, the mobility, which characterizes the linear response of the velocity to a constant 
external force /p(t) = /p as 

^^limM^, (21) 

is similarly determined through the relation /i = i?"'"(0). 

The bare friction coefficient 7, which appears in the model equation ^ can be determined 
from the high-frequency plateau of i?+(a;), because it satisfies 7 = lim^^^oo 

One may notice that the two plateaus are separated at an intermediate frequency, (indi- 
cated by an arrow in Figs.El(a)). This frequency corresponds to the decay rate of the probability 
distribution inside a potential well. Roughly speaking, this rate is the inverse of the time needed 
for a particle to diffuse around the bottom of the potential well. More precisely, this frequency 
coincides with the smallest absolute value of non-vanishing eigenvalues of the Fokker-Planck 
operator given in Eq. @. 

Note that these correspondences are not restricted to a system in equilibrium. The above- 
mentioned quantities can be obtained for nonequilibrium systems in the same manner. 



3.2 Nonequilibrium case: Violation of Fluctuation-Response Relation 

Next, let us examine what happens in a nonequilibrium case. We conducted the same calculation 
as above but with different parameter values: Q-iq{x) = ^oi{x) = 10 (7£^/T). In this case, the 
condition given in Eq. © is not satisfied and the canonical distribution is not the steady-state 
distribution. Then, from Eq. Q, the mean velocity does not vanish in general as mentioned ear- 
lier. Actually, if we employ the same potential profiles as above, Uq = sin(27rj;/£) + l/2 sin(47rx/^) 
and Ui = const., we numerically obtain v = 0.054 T/{^i). Furthermore, we find that the FRR 
[Eq. (fT^ ] is violated in this case as depicted in Fig. |21(b). In this figure, violation of the FRR is 
evident in particular for small frequencies. For the frequencies much larger than uj*, the FRR is 
retained even in a nonequilibrium steady state. This can be qualitatively explained as follows; 
For a timescale much shorter than l/w*, the particle does not know whether the potential is 
being switched because it does not travel a long distance within this short time range and the 
particle effectively feels no potential. Therefore, the statistical properties of the motion of the 
particle in the presence of the switching potential become same as those of a free diffusion of 
a particle for these short time scales, and the FRR holds. However, for a longer timescale, the 
particle starts to feel the potential as well as its switching. Then, the statistical property deviates 
from that of equilibrium system, and the FRR is violated. In the next Section, we consider the 
physical significance of violation of the FRR. 



4 Main Claim 

In this Section, we state our main claim. In the previous Section, it is found that the FRR is 
violated in a nonequilibrium steady state. Such a nonequilibrium steady state is maintained by 
the external energy input and energy dissipation into the environment. Then, by comparing the 
amounts of energy flows, J [see Eq. Q], with the degrees of the FRR violation in nonequilibrium 
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steady states, we are led to a simple relation between these quantities. This relation is expressed 
by the following equality 



This equality is proved in Sec. |SJ 

Interestingly, the right-hand side of Eq. (|22|) expresses the extent of the FRR violation. The 
first term in the braces, the square mean velocity, deviates from zero when the detailed-balance 
condition is violated in nonequilibrium states. The second term in the braces represents the 
extent of violation of Eq. H19|) over the frequency domain. This quantity is precisely equal to 
the area of the region surrounded by two curves C{uj) and 2k-QTR^{ijo) in a plot of them as 
functions of oj [see Fig. 121(b)] . Notice that the square-mean velocity can be interpreted as the 
zero-frequency component of the velocity correlation, (x(t)i;(0))g. Therefore, the right-hand side 
of Eq. (|22|) is interpreted as the sum of dissipated energy per unit time through each frequency 
of velocity fluctuations. In equilibrium, since the mean velocity vanishes and Eq. (|19() holds, 
the right-hand side of Eq. (|22j) vanishes. This is consistent with the fact that the rate of energy 
input and dissipation, J, vanishes in equilibrium. 

One of the outstanding features of Eq. (|22|) is that it is independent of the system details such 
as the profiles of the potentials and the transition rates. Indeed, this equality holds regardless of 
these model parameters. It is also independent of the number of the internal states. Furthermore, 
this equality can be proved for a wide variety of Langevin models other than the flashing ratchet 
model. It can be extended to time-dependent systems and systems which include multiple heat 
reservoirs. Therefore, this type of relation between the energy dissipation rate and the extent of 
the FRR violation is universal for a class of Langevin systems. See Ref. [20] for further details. 

We now numerically demonstrate Eq. (|22j) for the flashing ratchet model. In Fig. |31 the 
energy input rate J calculated as in Eq. Q and the right-hand side of Eq. are plotted 
independently as functions of the transition rate [The transition rates were taken uniform as 
VL = r2io(a;) = r2oi(3^)]- As clear seen in this figure, Eq. ()22|) holds for all the values of the 
transition rate O. 

5 Mathematical Proofs 

In this Section, mathematical proofs of the result described in the former Sections are presented. 

5.1 Path-integral representation of the flashing ratchet model 

As a preparation, in this Section, we develop a path-integral representation of the flashing ratchet 
model expressed as 




(22) 




(23) 



where 




(24) 



In this Section, the Boltzmann constant is normalized to unity. 
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Figure 3: Numerical verification of Eq. (|22|) for the flashing ratchet model. The solid curve 
represents the values of J calculated from Eq. @ . The filled squares are values of the right-hand 
side of Eq. H22|) . The horizontal axis represents the values of the transition rates fi. All the 
quantities are dimensionless under the normalization as £ = ^ = k^T = 1. 



5.1.1 Definitions and Notations 

Let us consider a standard description employed in quantum mechanics. We define the position 
and momentum operators as 

X = X (25) 
and denote eigenvectors of these operators as 



X\x) = x\x) (27) 
V\p) = p\p). (28) 

These eigenvectors are orthonormalized 

{x\x') = 5{x-x') (29) 
{p\p) = S{p-p'), (30) 

and regarded as complete sets in the state space Sj in the sense that 

;>oo 

dx\x){x\= / dp\p){p\=T, (31) 



where I is an identity operator. We can also confirm the following relations 

{x\p) = -i=e-^ (32) 



{p\x) = -^e-«^. (33) 
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Next, we consider a vector space, C^, and the following complete set in 



|0) ^ 

ID = 



(34) 
(35) 



We then consider the product space x Sj, whose element can be expressed as \(t,x) = \a) \x) 
with a G {0, 1}. \a, x) spans a complete set in x as 



dx(|0,x) (0,a;| + \l,x) (l,x|) =1. 



(36) 



Using the above notations, C{t) defined in Eq. 1)241) can be expressed as an operator acting on 
X ^ as follows 



at) = -v[u;,{A:)\o){o\ + u[{x)\i){i\]--Ut)v--v' 

'TO' 

-Oio(^) (|0) (0| - |1) (0|) + noi{X) (|0) (1| - |1) (1|) 



(37) 



5.1.2 Representation 

Next, we obtain a representation of C{t) in terms of the position. We first use the following 
identity 

C{t) = ffff (pil C{t) \p2) {x2\ ^e'("iP^-"2P2)dxidx2dpidp2. 



2tt 



(38) 



The following transformation of variables 



i(Xl + X2) 



X 

y = X1-X2 

V = \{Pl+P2) 

L 9 = P1-P2 



(39) 



yields 



where 



C{t) 



L(x,p, t)— e^^^dxdydp, 
27r 



L{x,p,t) = 
From Eq. H40() . we obtain 

{x\C{t) \x') 



P- 



£{t) 



dp 
— ( 
27r 



ip(x—x') 



L 



X + x' 



,p,'t 



(40) 
(41) 

(42) 



By substituting Eq. (|57|) into Eq. (|1T|) . it is calculated 



L{x,p,t) 



pU'o{x) - -U'^{x) 



|0) (0| + 



pU[{x) - -U'{{x) 



le 



|1) (l| --/p(t)p 

7 



- oio(x) (|o) (o| - |i) (o|) + noiix) (|o) (i| - |i) (i|) 

7 



(43) 
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5.1.3 Path-integral formulation 

Next, we integrate Eq. H23() from t = tj over a small interval At as 



Pi{x,tj+i) 



[I + AtC{tj)] 



(44) 



where tj+i = tj + At. Then, we insert Eq. H36|) with x = Xj+i and x = xj into the left and the 
right of Z^(tj+i, tj), respectively. Fm'ther integration and successive insertions of Eq. (|3()j) lead to 
the following expression of the probability distribution 

Par,{xN,tN) = W / dXj ^ {aj + l,Xj+i\U{tj + l,tj)\aj,Xj) Pa^,{xQ,tQ), (45) 

i=o •'-'^ <T,e{o,i} 

where terms of O(At^) have been omitted. From Eq. (|45|) . we obtain the transition probability 
of a path 

[cr, x] = {(o-o,Xo), (cJi,Xi), (crAr,XAr)} (46) 
which starts ((To,xo) at t = to as 



N-l 

Pe{{(TQ,XQ) \(T,X\) = Jl {aj+i,Xj+i\U{tj+i,tj)\aj,Xj) . 

3=0 

By use of Eqs. (|l2|) and (|43|) . we can calculate each factor entering Eq. (|47j) as 

(0,Xj+i|Z^(tj+i,tj) |0, Xj) = 



(47) 



7 



47rrAt 



exp 



At 
'4^ 



7 



Xj^l Xj 

At 



n 2 



+ C/^(x,)-e/p,, 



+ |^t^'(%)-^*^io(%)[,(48) 



where Xj = {xj + Xj+i)/2 and /pj = [/p(tj) + /p(tj_|_i)]/2. Similar calculations yield 
(l,Xj+i|Z^(tj+i,tj) |l,Xj) = 



7 



47rrAt 



exp 



At 
'4^ 



7 



At 



n 2 



+ C/{(x,)-e/p,^ 



+ |^t^"(%)-Atl7oi(x,-) 1,(49) 



and 



(1, Xj+i| Z//(tj+i, tj) |0, Xj) = Atr2io(xj)(5(xj-(.i — Xj), (50) 
(0,Xj+i|^^(tj+i,tj) |l,Xj) = I\tVlQi{xj)5{xj+i - Xj). (51) 

Using the above expressions, we can explicitly write down the transition probability. More 
concretely, let us suppose that the state a is switched at time intervals with j = ji,j2, ■■■,jn as 
Wjm+i) = (? o) I'^im) for m E {l,2,...,n}. Then, aj for jm < j < jm+i {rn E {0, l,..,n} and 
jo = —1, jn+i = N) has the same value represented by Sm- The transition probability of this 
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path is represented as 



ii V A-^TAf 



. 47rrAt 



X 

m=l 

/ ArT ll— T' _i2. 

J-J- V 47rrAt ' ^ ' 

j=jm + l 

where tm = -(Sm - !)• 

5.2 Proof of the fluctuation-dissipation theorem 

In this Section, we prove the fluctuation-dissipation theorem for the flashing ratchet model with 
aid of the results obtained in Sec. 15.11 

5.2.1 Local Detailed-Balance Condition 

First, we consider the time-reversed path of a path [u, x], denoted by 

[a,x] = {{aN,XN), (cTAr_i,XAr„i), {ao,xo)}. (53) 

Let fp{t) represent the time-reversed counterpart of fp{t) as fpj = fp^N-j-i- We then define the 
transition probability of a path [a, x] subject to fpj, which is represented by Ps{{(7o,xq) [a, x]), 
by replacing fpj in Eq. with fpj. With these preparations, we find 

Peiiao,Xo) [a,x]) ^ ^f}U,,^{xo)-(3U.^{xM)+PeY.^J^^ hA^j+i-^j) 
Pe{{aN,XN) [cr,x]) 

where f3 = T~^. If the condition in Eq. © is satisfied, Eq. ((51]) is reduced to 

Pe((o-o,Xo) [a,x]) ^ g/3C/<.o(xo)-/3C/<.j^(xjv)+/3eEf=V /pj(^i+i-^i). (55) 
Pe{{aN,XN) W,x]) 

This expression is a key property used in the proof of the fluctuation-dissipation theorem, and 
is termed the local detailed-balance condition j^HHH]- particular, by setting e = 0, Eq. 
can be transformed as 

P^^(xo)Po((^^o,xo) ^ [a,x]) = P^^{xN)Po{i^N,XN) ^ W,x]). (56) 

This relation is termed the detailed-balance condition with respect to the canonical distribution 
given in Eq. 0. From this expression, it can be concluded that the canonical distribution is 
the stationary distribution in the system when Eq. is satisfied. However, when Eq. (6) is not 
satisfied, there is no stationary distribution with which the detailed-balance condition holds. 
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5.2.2 Fluctuation-Dissipation Theorem 

We investigate the ensemble average of a quantity j4[(T, x], which is a functional of a path [o", x], 
when the initial distribution is the canonical distribution P^{x). The average is defined as 

N 



/oo 
dxj Yl Paoixo)Pe{i^o,xo)^[a,x])A[a,x]. (57) 

By using Eq. H55() and transforming variables from [a, x] to [a, x] , this expression can be rewritten 



as 

TV 



{A[a,x])^ = flTdxj Yl P^oi^o)Pe{{<To,xo)^[a,x])e-^'^^^oUA->.+i--^)A[a,x] 

-•r. J —OO ^ rr, n 



/N-l 



^/p,,(x,-+i-2;,)i[a,x] \ +0{e^), (58) 

where A[a, x] = A[a, x], and (• • •)^ denotes the average over processes subject to the time-reversed 
perturbation fp{t). 

We now suppose that N is an odd number and choose ^[fj, x] = (a;(^+i)/2 ~ X{N-i)/2)/^i- 
By setting e = in Eq. H58|) . we immediately obtain 



At 







Furthermore, we choose the perturbation such that fpj = fpj = fp^N-j-i- In this case, the two 
averages, (• • •)^ and (• • ■)-, coincide. Then, Eq. yields 

V (Af-3)/2 

=(TV+i)/2-X(^.,)/,\ ^ ^ ^ /3C^_/,,At + ./?^/„(^_,)/,At + 0(e^), (60) 



At 

j=0 



where 



is the velocity correlation function in the discretized form. The fact that Cj-k only depends on 
the difference j — /c is a consequence of the time-translational invariance of the system in a steady 
state including equilibrium. Cj-k converges to C{t) defined in Eq. (|12|) in the limit At ^ and 
N ^ oo with A^At fixed (hereafter, referred to as "the continuous-time limit"). 
The definition of the response function [Eq. 1)14(1 ] can be discretized as 

(7V-l)/2 

= Yl ^iv_i„./p,,At + 0(At), (62) 

j=0 

where Rj^t = R{tj — tk + At/2). In the first line in Eq. ()62|) . we have employed Ito's scheme in 
discretizing the right-hand side of Eq. (|14() . In what follows, we will use this scheme, although 
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the other schemes such as Stratonovich's scheme will lead to the same result if it is consistently 
used. 

Comparison of Eqs. (jSUj) and 1)62(1 yields 

r (3Cj for j > 1 
Rj = } (3Cq/2 for j = (63) 
otherwise 

Therefore, we finally obtain the desired result Eq. ()18|) in the continuous-time limit. 
Next, we introduce the discrete Fourier transform of a series {Aj} as 

(7V-l)/2 

^ AtA.e^™^ (64) 

j=~{N-l)/2 

which converges to the Fourier transform A{lo) with uo = 27ra/At in the continuous-time limit. 
Using this notation, we calculate 

(JV-l)/2 

= AtiJj-e^'^*"^' 
i=o 

{N-l)/2 

Y At/3Cje2™^' + At/?Co/2, (65) 

and 

-1 

R-^= Y Ai/?Cje2™^' + At/3Co/2, (66) 

j=_(jV-l)/2 

where the symmetry property Cj = C-j has been used. Combining these expressions, it is found 
that 

f3Ca = Ra + R^a- (67) 

Therefore, we obtain Eq. ()19|) in the continuous-time limit. 
5.3 Proof of the main claim 

In this Section, we prove our main claim Eq. (|22|) for the flashing ratchet model. First, we 
consider the steady-state average of a quantity A[a, x] = {xk+i — x^)/ At by replacing P^{x) with 
P|*(x) in Eq. ((57|) . In particular, the perturbation is chosen as fpj = 5ji/At. Next, Eq. (|57j) is 
rewritten in the following way 

C-^"^) = n n /dx.</dx.> E /dx./dw-(xo):^^^ 

^Pe((cro, xo) [a, x; ai = 0, cj/+i = 0]) 



(68) 



+Pe{{aQ,XQ) - 


^ [a, x; ai 


= 0, cr^+i 


= 1]) 


+Pe{{(JQ,XQ) - 


^ [o-, x; (7/ 


= l,cr«+l 


= 0]) 


+Pe{{ao,Xo) - 


■> [a, x; (7/ 


= l,cr«+l 


= 1])" 
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Then, by substituting Eq. into this expression and differentiating it with respect to e, it is 
calculated 



d_ / xk+i - Xk 
ds \ At 



l-l N 



\2lT 



J<=»]> 

- Xi 



At 



7- 



2-fT 



At 

xi+i - xi , 



Po((o-0, xq) [a, x; ui = 0, crz+i = 0]) 



At 



+ U[{xi) 



Po{{(yo,xo) [a,x]Oi = 1,0-^+1 = 1]) 



27T 



xi+i - XI , 
7 T7 + U„Xxi) 



At 



Xk+l - Xk 

At 



(69) 



Comparing this expression with Eq. (|62|) in which we set /pj = 6ji/At, we obtain the following 
relations 



v^ + jCk-i + (u'ixi 



Eq. (f7n|) can be reexpressed as 



Xk+l - Xk 

At 



Rk-i for I < k 
ioT l> k 



(70) 



Jj-k 

for k < j. It is then found 
Jo = - 



At 



7 (w^ + Cj_k - iTRj_k) 



(71) 



(f^a(i)(a:(t))oi(t)>Q = J, 



(72) 



in the continuous-time limit. In the same limit, the right-hand side of Eq. (|71j) converges as 

f ("^ r / \ 1 do; 

7 (t;2 + Co - 7r^o) ^ 7 <^ + / C(uj) - T (r{uj) + R{-ujy 



2tt ' 



(73) 



where Fourier's integration theorem 

R{-At/2) + RiAt/2) 
lim 

At->0 



^(w) + R{-uj) duj 



(74) 



2 7_oo 2 27r 

has been used. By combining these expressions, we finally obtain the desired result in Eq. (|22|) . 



6 Discussions and Conclusion 

6.1 Applications to experiments 

In this Section, we discuss the possibility and significance of applying our result to experimental 
studies on molecular motors. 

As explained in Sec.|3 the energy input to the system per unit time is obtained using Eq. (|22|) 
by determining the mean velocity, the power spectrum density, C(u;), and the susceptibility, 
R{uj), through measurements. This procedure requires no more detailed knowledge on the system 
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such as the potential profile. Therefore, if this method can be used in experimental studies of 
molecular motors, it will be a great advantage, since experimental determination of energy input 
rate is technically quite difficult. However, there are several restrictions in order to warrant the 
relevance of this method to a real molecular motor system. 

First, a motor system to be considered has to be reasonably well described with some kind of 
Langevin model including a Brownian motor model, since the present theory is constructed on 
the framework of a Langevin equation. This means that the position of the center of mass of the 
motor, x{t), has to be a considerably slow variable compared with the other degrees of freedom 
in the system. If there is another degree of freedom, y{t), that has a comparable timescale with 
that of x{t), we have to include y{t) in the system of Langevin equations to be considered. In this 
case, the quantity J given in Eq. (|22|) provides the information how much energy is dissipated 
through the degree of freedom x{t) per unit time. The total energy input will be the sum of all 
the dissipated energy through degrees of freedom x{t) and y{t). 

This condition might be satisfied in some cases. In particular, since an Fi motor has a single 
degree of freedom for rotation, it is plausible that this molecule can be described using a Langevin 
equation where the slow variable is a rotational angle of the rotational unit (7 subunit). A kind 
of kinesin that has only one head part, termed KIFIA [2^], is another example, since the motion 
of this molecule could be described with a Langevin equation where the slow variable is the 
position of the center of mass. Furthermore, a conventional kinesin that has two heads can also 
be modeled by a Langevin equation with a single variable, which is the position of the center of 
mass, because this molecule uses two heads alternatively and thus a single variable is enough to 
represent the state of the molecule. 

The most important point here is that the validity of these conjectures can be experimentally 
determined through examination of Eq. (|22j) . Because we can roughly estimate the free energy 
consumed by a single molecule per unit time by a macroscopic biochemical experiments, we 
can determine whether Eq. p2|) is satisfied in a molecule under consideration by adopting this 
macroscopic estimate of the energy input rate as J. If Eq. (|22|) is verified in experiments, 
this implies that our assignment of slow variables was correct. On the contrary, if it is denied 
experimentally, this implies that there are other slow variables and our knowledge on the system 
is insufficient. In this manner, experimental examination of Eq. H22|) can be used as a "check 
sum" to single out the slow variables that is involved in the energy transduction. 

Second, in order to correctly determine the right-hand side of Eq. (|22() . the time and spatial 
resolutions in experiments have to be sufficiently good, although very good time resolution is 
often difficult to be achieved. However, by noticing that the integrant of the second term on the 
right-hand side of Eq. (|22j) vanishes for high frequencies as seen in Fig. [2l (b), it is enough to 
determine C{uj) and R{uj) up to a finite frequency w* (the decay rate). From the investigation of 
the flashing ratchet model, it is conjectured that I'k joj* is close to the time needed for a motor 
to diffuse around a single potential well. Considering the current standard of the single-molecule 
techniques, this requirement seems to be within the scope. 

6.2 Example: KIFIA 

Here, we attempt to apply the present theory to KIFIA as an illustration. Unfortunately, the 
experimental data currently available in the literature of KIFIA is insufficient to determine C{ui) 
and RiyUj) in a sufficiently large frequency domain. Thus, we cannot verify the validity of Eq. ()22[) 
for this molecule at present. 

Instead, we here determine the decay rate lo* of KIFIA by assuming that the KIFIA can 
be modeled by the flashing ratchet. For this purpose, we employ an approximation for the 
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right-hand side of Eq. (|22l) . In the calculation of the FRR violation shown in Fig. [2 (b), it 
is found that the function C{uj) — 2k-QTR^{Lv) is well approximated by a Lorentzian, [C(0) — 
2k-BTR'^{dy\{oj*)'^ /[uj"^ + (tij*)^]. Thus, using this approximation in the second term on the right- 
hand side of Eq. (|22|) , we calculate 

/oo J 
C{u:)-2kBTR+{uj) ^(D-fikBT)uj*, (75) 
-oo L -I ^TT 

by use of the relations 2D = C(0) and fj, = R^{0). 

We then adopted the following values from Refs. [1311111111: v = 200 (nm/s), D = A x 10^ 
(nm^/s), = 1.3 X 10^ (nm/s-pN) in the presence of 2 mM ATP at T = 300 (K), and 7 = 1 x 10""^ 
(pN-s/nm). The value of 7 is determined as follows; when a KIFlA molecule is traped in a ADP- 
binding state, it freely diffuses on a microtubule with a diffusion coefficient -Dadp without energy 
consumption. We determined the value of 7 through Einstein's relation as 7 = /cbT^Z-Dadp by 
assuming that this is a free diffusion. 

As the value of J, we adopt the rate of free energy consumption through ATP hydrolysis 
reactions. In physiological condition, it is believed that a single hydrolysis reaction of ATP 
molecule produces the free energy of AG = 80 (pN-nm). A KIFlA molecule hydrolyses ATP 
molecules at the rate v = 1 x 10^ (s~^) in the presence of 2 mM ATP. Therefore, the energy 
consumption by a single motor per unit time is estimated as vAG = 8 x 10^ (pN-nm/s). 

We then determine the value of uj* through 

^ = ,iD-,k,Ty ^''^ 

Calculation using above mentioned experimental data yields lo* /27r = 4 x 10^ (s^^)- This value 
of the decay rate is not much different from the inverse of the cycle time of the ATP hydrolysis 
reactions and the dwell time of a KIFlA, it should be possible to determine the relevance of 
this value by use of the single- molecule experimental techniques. If the validity of the above 
obtained value of u* is experimentally verified, this implies that the flashing ratchet model is 
quantitatively relevant for the KIFlA. 



6.3 Concluding Remark 

In conclusion, a new equality relating the rates of energy input and dissipation to the extent 
of the FRR violation, which has been obtained by the authors recently, has been explained in 
detail. This relation enables us to investigate the energetics of a Langevin system on the basis 
of experimentally measurable quantities alone. With regard to the flashing ratchet model, the 
FRR and this new relation is numerically demonstrated and mathematically proved. It was also 
explained that this new relation can be used in experimental studies of molecular motors in order 
to single out the slow variables related to the energy transduction. An attempt to apply this 
theory to the KIFlA has been presented and the value of the decay rate has been predicted. 

We believe that to test this result for various types of molecular motors is important to 
elucidate the mechanism of their mechano-chemical energy transduction. Our method provides 
a novel approach in the study of molecular motors since it enables determination of energy flows 
at the single-molecule level. We also hope that by extending the present approach, we can 
establish a methodology to extract information of nonequilibrium systems including biological 
molecular machinery by external measurements. 

The authors acknowledge H. Noji, R. lino, K. Tabata, and K. Hayashi for discussions on 
several issues in this paper. The authors are also grateful to G. C. Paquette for his advices to 
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